Ray Tracing Simulations of Weak Lensing by Large-Scale Structure

We investigate weak lensing by large-scale structure using ray tracing through N-body simulations. Photon trajectories are followed through high resolution simulations of structure formation to make simulated maps of shear and convergence on the sky. Tests with varying numerical parameters are used to calibrate the accuracy of computed lensing statistics on angular scales from ∼ 1 arcminute to a few degrees. Various aspects of the weak lensing approximation are also tested. We show that the non-scalar component of the shear generated by the multiple deflections is small. For fields a few degrees on a side the shear power spectrum is almost entirely in the nonlinear regime and agrees well with nonlinear analytical predictions. Sampling fluctuations in power spectrum estimates are investigated by comparing several ray tracing realizations of a given model. For survey areas smaller than a degree on a side the main source of scatter is nonlinear coupling to modes larger than the survey. We develop a method which uses this effect to estimate Ωm from the scatter in power spectrum estimates for subregions of a larger survey. We show that the power spectrum can be measured accurately on scales corresponding to 1− 10 hMpc with realistic number densities of source galaxies with large intrinsic ellipticities. Non-Gaussian features in the one point distribution function of the weak lensing convergence (reconstructed from the shear) are also sensitive to Ωm. We suggest several techniques for estimating Ωm in the presence of noise and compare their